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ABSTRACT 

We present results of a semi- analytic model (SAM) that includes cold accretion and a porosity-based 
prescription for star formation. We can recover the puzzling observational results of low V/a seen in 
various massive disk or disk-like galaxies, if we allow 18 % of the accretion energy from cold flows to 
drive turbulence in gaseous disks at z = 2. The increase of gas mass through cold flows is by itself not 
sufRcient to increase the star formation rate sufficiently to recover the number density of M^, > 120 
Mq yr~^ galaxies in our model. In addition, it is necessary to increase the star formation efficiency. 
This can be achieved naturally in the porosity model, where star formation efficiency scales oc a, which 
scales as cloud velocity dispersion. As cold accretion is the main driver for gas velocity dispersion in 
our model, star formation efficiency parallels cold accretion rates, and allows fast conversion into stars. 
At 2: ~ 2, we find a space density 10~^ Mpc~^ in star-forming galaxies with M, > 120 Mq yr"-^, in 
better agreement than earlier estimates form SAMs. However, the fundamental relation between M* 
and is still offset from the observed relation, indicating the need for possibly more efficient star 
formation at high z perhaps associated with a role for AGN triggering. 
Subject headings: galaxies: evolution - galaxies: kinematics and dynamics - galaxies: spiral 



1. INTRODUCTION 

Multiple lines of evidence suggest that galaxy for- 
mation is taking place with a varying degree of ef- 
ficiency over cosmic time. The overall star forma- 
tion density i n the universe declines at z < 2 (e.g. 
iHopkind |2004| . and references therein), a feature seen 
also in the normalization of the rel ation between stel- 
lar mass and star formation ra t e (ICaputi et al.l 120061 : 
iDaddi et al.] l2007t iNoeske et all l2007f) . Further sup- 
port for the notion of rapid early formation of galax- 
ies at 2 > 2, especially massive ones, comes from high 
g/Fe ratios (iThomas et al.l2005l) , early metal enrichment 



(iMaiolino et al 



20081), grow th of super-massive black 

holes (jHasinger et al.' 2005' ) and the a ssembly rate of 
massive galaxies (Perez-Gonzalez et al.l [2008 ). The red- 
shift range 1.4 < z < 2.5 has been targeted by various 
surv eys combining broa dband color selectio n techniques 
(e^ lDaddi et al.l 120041 : iSteidel et all l2004 iDaddi et all 
l2007l (D07)) and recen tly by high-resolution Ha i ntegral 
field spectroscopy (|Forster Schreiber et al.l [2"006[ ). The 
latter reveals massive ~ IO^^Mq, rotating disk galaxies 
with star formation rates of ~ 140 Mq yr~^ and veloc - 
ity dispersions of ~ 60 km s~^ (e.g. iGenzel et al.ll2006l ). 
The origin of the high star formation rates and the large 
velocity dispersion in these galaxies is not presently un- 
derstood. 

On the theoretical side, detailed comparisons between 
observations at z = 2 and N-body-SPH simulations or 
semi-analytical models (SAM) show discrepancies. Al- 
though the overall slope of the M^—SFR is well recovered 
in the models, th e normalization is significantly smaller 
(DOT: lDavel[2008l see however Dekel et al. 2008). 

Theory and numerical simulations show that galax- 
ies accrete cold gas in two modes. At early times 
and in small dark matter halos, most of the gas 
never shock heats to the viria l temperature and is ac- 
creted on a free fall time scale (|Birnboim fc Dekel 120031 : 



Keres et al."2005; 'Dekel fc BirnboimI 120061 : lOcvirk et all 
2008:,Keres et al. 2009 ). The general condition for a sta- 
ble extended shock to be maintained is that the post- 
shock gas remains highly pressurized. In practice, this 
means that the loss of energy by radiative cooling needs 
to be at least compensated by the compression of post- 
shock infall material (Dekel fc Birnboim 2006 ). This 
condi tion leads to a criti cal halo mass oiMghock ~ 10^^-^ 
M0 (jOcvirk et al.l l2008f l for shock heating of gas and 
hence a transition from fast and efficient "cold accre- 
tion" onto the galaxy to a more moderate "hot accre- 
tion " mode set by the radiative coolin g time of the gas. 
In an earlier study iGenzel et al.l (|2008[ ) showed that the 
number density and accretion rates of dark matter halos 
suffice to provide enough high star forming galaxies, as- 
suming that the newly accreted dark matter provides the 
universal fraction of baryons, and that these baryons are 
all transformed into st ars. Alternat i vely, u sing a hydro- 
dynamical simulation, iDekel et al.l ()2009l ) showed that 
the gas accretion rates are comparable to the observed 
star formation rates thus requiring a very efficient trans- 
formation of accreted gas into stars. 

In this Letter, we aim at addressing the interlinked 
questions of high star formation rates and velocity dis- 
persion observed in galaxies at z = 2. Our work expands 
on the results previously presented bv lDekel et al.l ()2009l ) 
and Genzel et al. (2008). In particular, we model the ef- 
fect of cold accretion on driving turbulence in the inter- 
stellar medium (ISM) and employ a porosity description 
for the ISM and star formation within it. The latter hy- 
pothesis is able to provide high star formation efficiencies 
in the context of a physically motivated model. 

2. MODEL 

In the following, we briefly outline our fidu- 
cial semi-analytic m odel (jKhochfar fc BurkertI l2005l : 
iKhochfar fc Silkl l2006l . and references therein) and the 
changes we have introduced. We generate merger 
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trees for dark matter halos using a Monte-Carlo ap- 
proa ch based on the extended Press-Schechter formal- 
ism ([Somerville fc KolattI 119990 down to a mass resolu- 
tion of Mmin = 5 X 10^. Once a tree reaches Mmin, we 
start moving the tree forward in time including physical 
processes associated with the baryons within each dark 
matter halo that include gas cooling, star-formation, su- 
pernova feedback, reionization and merging of galax- 
ies on a dynamical friction time-scale. In this Letter, 
we will omit including detailed prescriptions for AGN- 
feedbac k fe.g. |Schawinski et al.|[2006.) or environmental 
effects (iKhochfar fc Ostrikei]l2noa[K 

In our modeling approach, we try to include the cold 
and hot accretion modes by first assigning the fraction 
fcoid of gas within the halo, that is in a cold accretion 
flow. If Mfiaio,g is the mass of gas within the halo, then 
we can readily calculate the amount of gas in the hot 
phase as (1 — fcoid)Mhaio,g- We use an approximate fit 
to t he fraction of cold a ccreted material as shown in Fig. 
3 of iK eres et al.' (2009) and the bottom part of Fig. 4 
in lOcvirk et al- (2 008). Please note that the results of 
iKeres et al.l ()2009l ) indicate a mild evolution with red- 
shift in this relation, which we for the sake of simplicity 
neglect in this study. Throughout this work we will use 
fcoid = -1.6 + 0.891og(M„,) - 0.06(log(M,„))2, with a 
lower and upper floor of and 1 for halos greater ^ 10^^'^ 
M0 and less ~ lO^^'^"* Mq, respectively. In a second 
step, we then calculate the cold and hot mode accre- 
tion onto the galaxy via Mcooi.c = fcoidMhaio,g/tdyn and 

Mcooi.h = O-H^ - fcoid)iMhotrcooi)/itdynRvir), respec- 
tively. Part of the underlying assumption is that cold 
accretion flows will reach the galaxy on a halo dynamical 
time tdyrn and that the hot accretion rate will be gov- 
erned by the radiatively coo li ng tim e using recipes laid 
out as in e.g. ISpringel et al.l ()2001l ). We find in general 
that cold accretion clearly dominates over hot accretion 
in agreement with results of numer ical simulations (e.g. 
iDekel et al.]l2009t IKeres et al.ll200 9l). The accretion rate 
of gas onto the disk of the central galaxy is then approxi- 
mated as Mcoid — Mcooi.c + Mcooi.h- This choice of accre- 
tion rate is the simplest possibility and shows results in 
agreement w ith those from hyd rodynamical simulations 
presented bvl Dekel et aT (2009). 

Star formation is modeled following the Schmidt- 
Kennicutt-law Af* = aMcoid/tdyn,disc (jKennicuttj 1 1 998f) 
with a porosity-based pres cripti on for the multiphase 
structure of the ISM (Silk 1200 ll ). Momentum balance 
yields a star formation efficiency a « 0.02((Tg/10 km 
s-^) (vcooi/400 km s-i) (msw/150 Mq) (lO^i erg/EsN) 
([Siig|200^ . The ag dependence arises because of the 
enhancement in star formation efficiency via momentum 
input from supernova remnants with increasing ambi- 
ent turbulent pressure. Note that gas accretion also 
drives gas turbulence, as in our model. Momentum input 
is only effective once the remnant expansion speed has 
slowed to below Vcooi ■ Here ag is the gas velocity disper- 
sion, msN is the IMF-dependent mass in stars formed 
for each type II supernovae, Esn the energy released 
by each supernova and Vcooi the shell velocity at which 
strong radiative energy losses set in. We set Vcooi — 400 
km s~^ and use a Chabrier IMF with msN = 150 Mq 
throughout this paper. To calculate the gas velocity dis- 
persion in the disk, we assume equilibrium between the 



rate at which energy is dissipated in the disk and energy 
added due to kinetic energy of infalling gas and super- 
novae. Including the cr dependence of the star forma- 
tion rate predicted by the porosity model, this leads to 

<j ^ (3 + {l3^ + vMcool,ctdyn,d^sc/McoldVY^^■ The free 
parameter 77 allows to adjust for our ignorance with re- 
spect to the actual fraction of the kinetic energy being 
used to drive turbulence and the relation between the 
dissipation time-scale and dynamical time of the disk, 
r/ is set to match the observed relation between V and 
V/a of galaxies at z ~ 2. The constant j3 = 10/2 km 
s~^ is set to ensure that in the limiting case of no ac- 
cretion, the local observed velocity dispersion in disks is 
recovered. Simulations show that cold accretion comes 
in along filaments on almost r adially orbits before hitting 
the galactic disk (jOcvirk et al..,2008. ). thus being able to 
contribute a large fraction of the gas kinetic energy to 
the heating of the gaseous disk, in contrast to hot-mode 
accreted material that settles more gradually onto the 
disk. 

In this study, we simulate a volume of 10^ Mpc^, and 
use the following set of cosmological parameters derived 
from a combination of the 5-year WMAP data with Type 
la supern ovae and measuremen ts of baryon acoustic os- 
cillations (jKomatsu et al.ll2008l ): = 0.28, = 0.72, 
nb/no = 0.16, cTg = 0.8 and h = 0.7. 

3. LOW V/a DISK GALAXIES 

We start by investigating the increase in gas velocity 
dispersion in gaseous disks due to cold accretion. Fig. 
1 shows the correlation between V/a and V. In gen- 
eral only 18% {rj = 0.18) of the accretion energy is suf- 
ficient to drive tu rbulence that agrees with the range 
of ob served values (|Genzel et al.ll2006l 120081: IStark et all 
|2008| ) . It is interesting to note that this value is in agree- 
ment with recent theoretical predictions by Dckcl et al] 
(|2009t) . Cold accretion-d riven turbulence in galaxies liv- 
ing in small dark matter halos with Vc < 100 km s~^ 
results in velocity dispersions similar to or larger than 
the rotational velocity, mainly because of the high ac- 
cretion rates and the low gas masses in the disks to 
absorb the associated accretion energy. These galaxies 
are dynamically hot system that are pressure-supported 
against gravitational collapse in contrast to classical cold 
disks. More massive galaxies with larger V typically live 
in massive halos which only have small fractions of cold 
accretion via filaments and massive gaseous disks that 
are not easily heated. The steep transition from low to 
high V/ a reflects the steep decline of cold accretion flows 
as a function of dark halo mass that we have adopted 
based on simulations. The scatter seen in the V/a — V 
plot is largely due to the varying accretion histories of 
individual halos. 

4. STAR FORMATION RATES AND MASS FUNCTION 

The inclusion of explicit cold accretion in the SAM 
increases the overall accretion rate onto model galaxies. 
This might suggest that the star formation rates will in- 
crease as well. However, using a constant star forma- 
tion efficiency normalized to the locally observed value 
of a = 0.02 the number density of star forming galaxies 
with Af, > 120 M0 yr~^ falls short by almost an or- 
der of ma gnitude in a greement with previous reports by 
D07 and (jPavel [20081 ) . FueHng of the galaxies with gas 
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Fig. 1. — The relation between V/cr and V for modeled galaxies 
(small filled circles). The observat ions are from lGenzel et a Tl ||2006|) 
(large filled circle), Genzel et al.l (2008) (filled stars). l^tarkeT^Ll 
|l008) (filled triangle) and Cresci (2008, ApJ submitted) (filled 
squares). We here assume that 18% of the accretion energy from 
cold fiows is used to drive turbulence in the disk. Note that mergers 
are excluded as well as satellite galaxies that by construction do 
not accrete cold gas anymore. 

is not the limiting factor for high star formation rates. 
Rather, it is the conversion of gas into stars over a short 
time-scale. The porosity-based model for star formation 
predicts that star formation efhciency scales as a oc ct, 
thus providing a natural way to increase the star for- 
mation rate with redshift in parallel with the cold gas 
accretion rate. 

In Fig. 2, we present the SFR-Af, relation for central 
galaxies at 2 = 2. The model galaxies cover most of 
the observed range. However, we find that the slope is 
slightly tilted. Mostly massive model galaxies have too 
low star formation rates. Inspection of the D07 data 
reveals that the statistics at the massive end is rather 
poor and that a possible flattening at high masses as seen 
in our models cannot be ruled out. It is interesting to 
note that the modeled SFR— A/* relation at 2 = 2.5 lies 
perfectly on the region covered by the observed galaxies. 
The number density of star-forming galaxies with > 
120 Mq yr-i at z = 2 is 10"'' Mpc"^ and within a factor 
two of the observed density of 1.6 x 10~^ in D07. Again 
our estimates for z = 2.5 are in excellent agreement with 
the results of D07 at z = 2. The number density of star- 
forming galaxies increases by a factor tw o to z = 3.5 and 
then starts to drop again (see Fig. 3). D ekel et al.l (120091 ) 
find in their simulation a similar increase in the number 
density of galaxies accreting gas with 120 yr~^ going 
from z = 2toz = 3 — 4. However, their overall number 
densities are larger by a factor ^ 2 or more. Also do 
they not find a peak around z ~ 3 as we do, but a steady 
decline in the number density of systems accreting with 
high rates. 

As we are not generating light cones to compare with 
the observations, but use redshift slices, we need to ask 
the question whether the ULIRG sub-sample of D07 is 
biased toward z > 2 within the overall sample which 
spans the range 1.4 < z < 2.5. Inspecting their Fig. 
6, we find that galaxies with L/^ > 10^^ Lq span a 
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Fig. 2. — Relation between star formation rates and stellar mass 
for central galaxies. The open squares show the D07 data. The 
pattern in the model data is due to the binning of halo masses 
along the halo mass function at the beginning of the merging tree 
at z = 2. 

wide range in redshift between z = 1.5 and z = 3 with 
brighter galaxies preferably being at z > 2, suggesting 
that we are underestimating number densities because we 
neglect the z > 2 contribution. Taking this into account, 
the difference with D07 should be even less than a factor 
of two. 

Finally, we turn to the overall mass function of galax- 
ies at z = 2, shown in Fig. 4. One might expect that 
due to the enhancement of early star formation in the re- 
vised model, basic properties like the mass function will 
change significantly. However, good agreement is found 
between the modeled mass function and those base d on 
Spitzer observatio ns bv | Perez-Gonzalez et al.l (|2005 ) and 
a recent study bv iMarchesini et all (120081) . The GOODS- 
MUSIC results of lFontana et"an ([20o¥ are a bit lower, 
but agree within their error margin. 

5. CONCLUSIONS 

We have investigated a possible scenario for the forma- 
tion of galaxies at high redshifts that is driven mainly by 
cold accretion flows. Massive galaxies at z = 2 live in ha- 
los that are at the intersections of fllaments that form the 
cosmic web. As shown in numerical simulations, these fil- 
aments provide a steady stream of dense cold gas that 
can reach the galaxy on a free-fall time scale and fuel 
ongoing star formation. Not only do these streams pro- 
vide cold gas, but they can also drive turbulence in the 
gaseous disk by depositing some fraction of the accretion 
energy into the ISM. We find that between 10% - 20% 
of this accretion energy is sufficient to increase the gas 
velocity dispersion to observed values. In particular we 
choose T] = 0.18, which still leaves accretion driven tur- 
bulence dominant with respect to supernovae driven one. 
A natural consequence of cold accretion driving turbu- 
lence is that it becomes less important in massive halos 
> lO"'^^'^ Mq and at low redshifts, where the fraction of 
cold mode accretion drops and V ja increases. Low red- 
shift observations find velocity dispersions of 10 km s"^, 
much lower than those at z = 2. This can be understood 
as a natural consequence of the steady demise of cold 
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Fig. 3. — Number density of star-forming galaxies with SFR 
> 60 M0 yr-i (stars), > 120 M0 yr"! (circles) and > 300 M0 
yr~^ (triangles) . The square is the data of D07 for galaxies with 
SFR > 120 Mq yr-i. 
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Fig. 4. — Galaxy mass functions at 2 = 2. Filled circles show 
results including cold accretion and the porosity based model for 
star formation. 



accretion with redshift in halos that host disk galaxies 
today. The host halos of these galaxies live no longer 
at the intersection of filaments, but in the filaments it- 
self. Those galaxies living in halos above the critical dark 
matter mass for stable shocks (Dckc l fc BirnboimI [20061 ) 
will no longer accrete from cold flows, but from hot halo 
gas, if at all, that radiatively cools and the main driver 
of turbulence will be supernovae. Such a picture does fit 
in nicely with observation of low redshift disk galaxies, 
which have much larger V/a than their high-z counter 
parts l|Kassin et al.ll2007[ ). 

Only increasing the amount of available fuel due to 
cold accretion is not enough to increase the number 
density of galaxies with large star formation rates. It is 
necessary to also increase the star formation efficiency 
at the same time. We tested an analytic model for a 
multiphase ISM that is based on a porosity prescription, 
that has been shown to capture some of the relevant 
physical mechanism related to star format ion seen in 
high-resolution simulations ijSlvz et al.ll2005[ ). The star 
formation efficiency in this model is a oc a, resulting 
in an increasing efficiency that parallels the gas supply 
by cold accretion and allows to convert it into stars 
on a very short time-scale. With this model we find 
improved agreement with the observed number density 
of star-forming galaxies. However our predicted star 
formation rates are lower than in the observed z ~ 2 
sample especially for the more massive galaxies. We 
have two comments on this point. Firstly, the data 
are most likely biased to high star formation rates. 
Secondly, galaxies with high star-forming rates at high 
redshift commonly show weak AGN activity (see for 
example a recent s tudy of a sample of z ~ 1.9 ULIRGs 
(jHuang et al.ll2009( l ). This most hkely is the tip of the 
AGN iceberg because of possible obscured AGN. AGN 
outflows could trigger as well as play the more accepted 
role of quenching star formation at high redshift. It 
remains to be investigated how well the agreement at 
different redshifts is, and we plan on addressing this in 
a future paper. 

We would like to thank the referee Avishai Dekel 
for his very helpful comments that helped improve the 
manuscript, Andreas Burkert and the SINS team for use- 
ful discussions, and Emanuele Daddi as well as Giovanni 
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